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A wide range of quasi-one-dimensional materials, consisting of weakly coupled chains, un- 
dergo three-dimensional phase transitions that can be described by a complex order parameter. 
A Ginzburg-Landau theory is derived for such a transition. It is shown that intrachain fluctua- 
tions in the order parameter play a crucial role and must be treated exactly. The effect of these 
fluctuations is determined by a single dimensionless parameter. The three-dimensional transition 
temperature, the associated specific heat jump, coherence lengths, and width of the critical re- 
gion, are computed assuming that the single chain Ginzburg-Landau coefficients are independent of 
temperature. The width of the critical region, estimated from the Ginzburg criterion, is virtually 
• parameter independent, being about 5-8 per cent of the transition temperature. 
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I. INTRODUCTION 
A. Motivation 

> 

OO 

\^ • Crystals consisting of linear molecular chains can have highly anisotropic electronic properties. Many of these 
' quasi-one-dimensional materials undergo finite temperature phase transitions. Detailed experime|ntal studies in the 
1970's of the charge-density-wave (CDW) transition in materials such as KCP and |TIF-TCNQ|J stimulated work 
' on the theoretical description of phase transitions in quasi-one-dimensional materialscTLl. Unfortunately, due to the 
limited quality of samples it was not possible to make a quantitative comparison of experiment with theory. In 
the past decade,p.a whole new range of quasi-one-dimensional materials, many in high-purity single crystals, have 
been synthesizejH. Examples of new materials are the CDW material K0.3M0O3 (blue bronze), the Bechgaacji salt 
(TMTSF)2PF60, that undergoes a spin-density-wave (SDW) transition, j-thin wires of superconducting leadlJ, and 
an inorganic compound CuGeOs that undergoes a spin-Peierls transitionlld. Furthermore, high-quality experimental 
datalljEJ on these crystals opens the possibility of making a quantitative comparison of experiment with theory. For 
f~\ ^ example, recently anomalies in the specific heat. Young's modulus, thermal expansion, and magnetic susceptibility 
O ■ close to the three-dimensional CDW transition in K0.3M0O3 were all precisely measured on the same single crystalEj. 
^ ' Fluctuations in the order parameter are extremely important in these materials. In a strictly one-dimensional 
^ , system with short-range interactions there a re M n t) phase transitions at finite temperature because fiuctuations in 
the order parameter destroy long range orderfijllj. Consequently in a real quasi-one-dimensional material a finite- 
^ . , temperature phase transition only occurs as a result of the weak interchain interactions. Hence, a complete theory 
must treat fluctuations along the chain carefully and also include the interchain interactions. On the other hand, 
in most materials the three-dimensional transition temperature, T^d, is clearly defined and the width of the critical 
region is only a few percent of T^jj. Hence, a three-dimensional Ginzburg-Landau theory should accurately describe 
the transition, except in a narrow temperature range. 

The purpose of the present paper is to provide a complete derivation of a Ginzburg-Landau theory describing 
the three-dimensional ordering transition of a quasi-one-dimensional system with a complex (i.e., two-component) 
order parameter. Such an order parameter describes superconductors, charge- and spin-density waves (which are 
incommensurate with the lattice), and a spin-Peierls system in a large magnetic field. Partic.ular attention is given 
to examining the validity of the-, assumptions and approximations made in previous worktTLl. This work is the 
continuation of a larger progranJllTEj of examining the effect that lattice fiuctuations have on electronic properties of 
CDW compounds. Recently, the role of ordet,parameter fiuctuations in the blue bronze and the Bechgaard salts was 
considered by Castella, Baeriswyl, and Makic3. Schulz and Bourbonnais-have considered quantum fluctuations in the 
phase of the order parameter in quasi-one-dimensional superconductor£il. 
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B. Overview 



The Ginzburg-Landau free energy functional Fi[(j3\ for a single chain with a complex order parameter 4){z), where 
z is the co-ordinate along the chain, is 



dz 



dz 



(1) 



In this paper the coefficients a, &, and c will be treatetLas phenomological parameters. For a specific system these 
coefficients can be calculated from microscopic theoryuc3. , ■, , 

Due to fluctuations in the order parameter this system cannot develop long-range order at finite temperaturetSllj. 
Furthermore, the values for the coefficients a, 6, and c given by the simplest microscopic theories predict that tlae 
fluctuations are important over a temperature range comparable to the single chain mean-field transition temperatureQ. 
To describe a finite-temperature phase transition, a set of weakly interacting chains are considered. If 4>i{z) is the 
order parameter on the i-th chain the free energy functional for the system is 



F\. 



.(z)]=^Fi[0,(z)] 



J,:, 



i,3 



dz'R.e[(j)i{z)* <j)j (z)] 



(2) 



where Jy describes the interchain interactions. In most of this paper it will be assumed that the interchain interaction 
Jij is non-zero only for nearest neighbour chains and that its value is Jx/"^ and Jj,/4 in the x and y directions 
respectively. A mean-field treatment of this functional will not give accurate results due to the large intrachain 
fluctuations (see Section III). This problem is solved by integrating out these fluctuations to derive a new Ginzburg- 
Landau functional with renormalizcd coefficients. 
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where Ux and ay are the lattice constants perpendicular to the chains. The new order parameter z), is 

proportional to the average of 4>i{z) over neighbouring chains (see equation (^2|)). The three-dimensional mean-field 
temperature T^d is defined as the temperature at which the the coefficient A{T) changes sign. Close to T^d 



\T3D 



(4) 



The goal of the present paper is to derive the functional (^ and find expressions for the transition temperature 
T^D and the coefficients A' , B, Cx, Cy, and Cz in terms of the interchain interaction and the coefficients a, b, and 
c of a single chain. 

The coefficients determine measurable quantities associated with the transition such as the specific heat jump, 
coherence lengths and width of the critical region. The specific heat jump per chain at the transition is given by 



l\2 



AC 



2BT3D ■ 



(5) 



The specific heat jump per unit volume is ACsd = AC / {axtty) . In the Gaussian approximation the correlation length 
parallel to the chains is given by 



-1/2 



where ^oz is the longitudinal coherence length defined as 



(6) 



(7) 



The correlation and coheMice lengths perpendicular to the chains are given by similar expressions. 
The Ginzburg criterionO gives a rough estimate of the width, AT^d, of the 3D critical region. 
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This paper establishes the striking result that most of the physics is determined by a single dimensionless parameter 



K=-^. 9 

lal'^c 

A brief argument is now given to show that k is a measure of the fluctuations along a single chain. The rms 
fluctuation in the single chain order parameter <| </> p>, calculated in the Gaussian approximation is T/2(c|a|)^/^. 
The magnitude of the order parameter, calculated in the mean-field approximation, (po, is given by 0q = |a|/26. 
Hence, k = 2{<\ (jj |2> /^g^^ 

so is a measure of the importance of fluctuations. The strength of the interchain 
interactions determines the value of k at the three-dimensiptiml transition temperature. 

The general approach that has been taken previouslyauo when deriving a Ginzburg-Landau functional for the 
three-dimensional transition is to treat the interchain interactions in the mean-field approximation and then solve 
the res ultin g one-dimensional problem treating the fluctuations along the chain in the lowest-level approximation (see 



Section A 2 for a definition) . In this paper the interchain interactions are also treated in the mean-field approximation 
but the fluctuations along the chain are treated exactly. This paper is confined to a mean-field analysis of the functional 
(^. In Section IV it is shown that this is justified except in a narrow temperature range very close to T^jj. However, 



it should be pointed out tiiat one could perform a sophisticated renormalization group analysis for the functional (|^) 
such as that due to Chenc3|^d which has been recently used in the analysis of specific heat measurements on the 
CDW compound K0.3M0O3E3. The free energy functional that derived here could be used as the input to such an 
analysis. 

The outline of the paper is as follows. Table ^contains a summary of the symbols for the important quantities for a 
single chain and for a three-dimensional system of weakly interacting chains. Section II describes an exact treatment 



of the fluctuations in a strictly one-dimensional system. Section [II contains a complete derivation of the coefficients 
in the Ginzburg-Landau functional (|^). In Section [V it is assumed that the single chain coefficients a, b, and c are 



temperature independent. The three-dimensional transition temperature is calculated as a function of the interchain 
coupfing (Figure Q). The values derived for the coefficients A', B, Cx, Cy, and Cz are used to calculate the specific 
heat jump, coherence lengths, and width of the critical region at the three-dimensional transition (Figures ^ and 
The striking result is established that the width of the critical region is fairly independent of any parameters, 
At3D 0.05 - 0.08. Sections g and |lH can be omitted by readers not interested in technical details. The Hartree, 



Hartree-Fock, and lowest-level approximations for the fluctuations that have been given in previous work are described 
in the Appendix. It is shown that the lowest-level approximation is only reliable for k < 10~^ (Figures i| and In 
Table |l| very rough estimates of k for various CDW materials give k > 10~^, showing that they are well outside the 
regime of validity of the lowest-level approximation. 

II. ONE-DIMENSIONAL GINZBURG-LANDAU THEORY 

The mean-field treatment of one-dimensional Ginzburg-Landau theory is described before an exact treatment of 
the order-parameter fluctuations is given. For completeness the Appendix to the paper describes various approximate 
treatments of fluctuations that have been given in earlier work. 

A. Mean-field treatment 

Long-wavelength fluctuations associated with the gradient term in (|^) are neglected. One simply minimizes the 
potential 

y(0) - a I p + 6 U r (10) 

The single chain mean-field transition temperature Tq is defined by the temperature at which the coefficient a(T) 
vanishes. Close to Tq it is convenient to write 

a{T)^{t-l)a' (11) 

where t = T /Tq. In the mean-field approximation there is a phase transition at To and the magnitude of the order 
parameter for T < Tqi '/'O; is 

'^o - (12) 
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The free energy per unit length of the ordered state, relative to the disordered state is smaller by Vq, given by 



Vo ^ -V[M = ^- (13) 



The jump in the specific heat at the transition is 



AC\d = (14) 
2 6 To ^ ^ 



An important length scale is the coherence length ^o, defined by 

The one-dimensional Ginzburg criterionil provides an estimate of the temperature range, ATio, over which critical 
fluctuations are important. 



There are several serious problems with a mean- field treatment. First, it predicts a finite temperature phase 
transition. Second, the simplest microscopic models give Ati^i ^ fl3, suggesting that the neglected fluctuations play 
an important role over a very broad temperature range. An exact treatment of these fluctuations is now given. 



B. Exact solution 



The goal is to evaluate the functional integral for the partition function Z of a system with the Ginzburg-Landau 
free energy (m). 



Z 



[d0(z)]exp(-Fi[0(z)]/r) 



It is convenient to rescale the order parameter 0(2:) —f (j){z)/ipo and write the functional as 



dz 



-2 



2c 



dz 



(17) 



(18) 



assuming that T < Tq. Scalapino, Sears, and FerrelO showed that a transfer matrix method could be used to reduce 
the problem to that of diagonalizing a transfer matrix Hamiltonian Htr- This involves a cylindrically symmetric 
Schrodinger-like equation with a single complex degree of freedom s — pe^"^ , 
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dp^ p dp p^ dip'^ 



-2p' 



*m(/0, <^) = Xm^mip, ■ 



(19) 



The eigenvalues Am only depend on the dimensionless parameter k. 

As discussed ia the introduction k is a measure of the size of the fluctuations in the order parameter, and defined 
by equation (^c3. k is central to this paper. It is shown in Section IV that the value of k at the three-dimensional 
transition temperature determines physically important quantities such as the specific heat jump. For the case of a 
real (i.e., one-component) order parameter the importance of k was emphasized previously by Bishop and KrumhanslEl 
(k corresponds to in their paper) and Dieterichtl {k corresponds to 4/q!^ in his paper). 

The partition function ( p7|) for a system of length L is 



Z = Y,cxp{-LVoXrn/T) 



(20) 



Consequently, in the thermodynamic limit, the free energy per unit length equals VqXq. The order parameter corre- 
lation function is 
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(mm*) = <^oE l< *™I^I^O >|' exp (-\^VoiXrn - Ao) 

which at large distances becomes 

{<Piz)m*) = 00 l< *ikl*o >P exp (^-M 



(21) 



(22) 



where is the correlation length given by the separation of the lowest two eigenvalues of the transfer matrix Hamil- 
tonian 



14)(Ai-Ao) 



The Fourier transform of the correlation function defines the static hnear susceptibihty 



For small q it is of the form 



where 



xiiq) = 



XI (0) 



Xi 



202^ |<^'„|s|*o>P 4 



Am — Aq 



a 



(23) 



(24) 



(25) 



(26) 



and a dimensionless function, /(k), that only depends on the fluctuation parameter k has been introduced. 
Due to the cylindrical symmetry of the potential V{s) the eigenfunctions can be written in the form 

«'„,,(p, ip) = ^^^^ exp(i£(^) £ = 0, ±1, ±2, ... 

Equation ( p^ reduces to the one-dimensional Schrodinger type equation 



k{£' 



2_ 1] 



dp2 



+ p^ 



UnAp) = XrijUn,l{p)- 



(27) 



(28) 



In order for the system wave function to be finite at the origin the radial wave function satisfies the boundary condition 

w„,HO) = 0. (29) 

Normalization of the wave function requires that 



I UnAp) ? dp=\ 



(30) 



The integrals over the angular variable ip in the matrix elements in the susceptibility (^6|) can be performed explicitly. 
The function /(k) in the susceptibihty ([2^) then reduces to 



/(-) = E 



|< UnAp\unfi >| 

An,l — Ao,o 



(31) 



This expression has been evaluated numerically by solving the eigenvalue equation ( |28| ) numericallyEj. The results 
are shown in Figure 0. To check the numerical results use was made of the sum rule 



^ '^V^ri,! - Ao,o) |< Uns\p\uO,Q >P= 1 



(32) 
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which can be derived using standard argument^. 

It is useful to consider the effect of an external potential $ on the system. Let A($) be the energy of the lowe 
eigenstate of the transfer operator Ht,- + s<i>* + s*<i>, where Htr is given by ([Tg|). Second-order perturbation theory 
shows that the static linear susceptibility xi(0) is given by 



x.(0)- *»^^^'*> 



The third-order nonlinear susceptibility is 

4 



X3 



(33) 



(34) 



" $=0 

Expanding the Brillouin-Wigner perturbation theory expression for A($)il to fourth-order in |<I>| gives 
^'~Kf (A™-Ao)(A„-Ao)(Ap-Ao) 

^^^^ (A™-Ao)(A„-Ao)(Ap-Ao) 

'2^1 \<^^\s\^o>\^ _ 16b 

- W''^'^ ? (A,„ - Ao)^ - -^'^'^ ^''^ 

As for the linear susceptibility the above expression can be simplified by using the factorization of the wavefunction 
( p7| ) and perfor ming the integrals over the angular variable Lp in the matrix elements. The dimensionless function 
gin), defined in (p5[), only depends on the fluctuation parameter k and is given by 



< U0fi\p\u,n,l X UmMp\u7i,0 >< Un,o\p\Up^i >< Up,l|jO|Mo,0 > 

(Am,i — Ao,o)(A„,o — Ao,o)(Ap,i — Ao,o) 

< Uofi\p\Um,l >< Uma|p|u«.2 >< UnMPhpA >< "p.iIpI^O.O > 



(A„i,i — Ao,o)(An,2 — Ao,o)(Apa — Ao,o) 



m,n,p 

"4/(.)E '^a"''''r^^'' - (3^) 

^ (A„,i - Ao.o) 

The functions /(k) and g{K) will appear in the expressions for the coefficients of the Ginzburg-Landau functional for 
the three-dimensional transition. 

In summary, the free energy per unit length in the presence of an external field $, -Fi[<I'], is, to fourth-order in |<I>| , 

= yoA(O) - XI (0) I $ P -X3 I $ 1^ (37) 
If the external field is slowly varying in space ^{z) then the quadratic term is replaced with 

Hz)* I dz'xiiz - z'Mz') = ^zT I ^ XiiqmiV 

dq 
2^ 



Hz)* I XI (0) {l-q^ei)^^^ 



= Xi(0)$(z)*^$(z)+e?^^) (38) 
where the long-wavelength expression (E3) has been used. 



III. DERIVATION OF A THREE-DIMENSIONAL GINZBURG-LANDAU FUNCTIONAL 

In this section a new three-dimensional Ginzburg-Landau free energy functional is derived by integrating out the 
order parameter fluctuations along the chain exactly and treating the interchain interactions in the mean-field approx- 
imation. Since this treatment (over)emphasizes the one-dimensional fluctuations its regime of validity is considered 
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briefly in the Conclusions section. Before proceeding with the derivation a rough estimate is made of how small the 
interchain interactions must be in order for a mean-field treatment of the anisotropie, functional (^) to be invalid. 
Taking the continuum limit of (||) and using the three-dimensional Ginzburg criteriorH gives 

M,D^\{M,nf^. (39) 

Hence, if Atuj ^ 1, as predicted by the simplest microscopic theoriesB, and if ^ Jy < 0.3a' (the regime considered 
in this paper) then At^jj > 1 and a mean-field analysis will give poor results. 

The treatment siijen here is similar to the derivation of the Ginzburg-Landau functional for the Ising model by 
Negele and OrlandEa. The partition function for a set of chains is 

Z = Ju^[dMz)]expi-F[Mz)]/T) (40) 

where the free energy functional F[<j)i{z)] is given by (|^). The interaction term is replaced with the following Gaussian 
integral 



I ^ X] '^y / Rc[(/),(z)*(/)j(2;)] 
= A /"n,[d$,(z)]exp -T^J^i / Re[$,(z)*$,(z)]+^ / 



dz Re[0,(z)*$,(z)] (41) 



where A is a normalization constant. The new field ^i(z) will be used below as the order parameter to describe the 
three-dimensional transition. It follows from (^0|) and ( p] ) that 

j 

Hence, the new order parameter is the average of the single chain order parameter <j)j{z) over neighbouring chains. 
The new expression for the partition function has the desirable feature that there are no interchain interactions 
involving the order parameter (t>i{z). Consequently, the results of Section |II B| can be used to integrate exactly over 
the fluctuations in (j)i{z) along each chain. The partition function is now a functional integral over the new order 
parameter $i(z) which has a free energy functional i^[<I>i(z)], 

Z = j n,[d$,(z)] exp (-F[$,(z)]/t) (43) 

where 

F[<I>.(z)]=^i?i[<i>,(z)]+T^J^iy' dzRe[$,(z)*$,(z)] (44) 

and Fi[$(z)] is the free energy of a single chain in the presence of an external field <i>(z), 

exp (-Fi[$(z)]/r) = J [d0(z)]exp J^-^iMfll +^ J dzRe[0(z)*<I)(z)]^ . (45) 

Up to this point the analysis is exact. A mean-field analysis of the functional integral (^) will now be given. In the 
next Section it is shown that this is justified except in a narrow temperature range very close to T^o- However, it 
should be stressed that one could perform a sophisticated renormalization group analysis of (^3|). The free energy 
functional that derived here could be used as the input to such an analysis. 

Close to the three-dimensional transition temperature the new order parameter ^i{z) will be small and slowly 
varying in space. Hence Fi[$(z)] can be evaluated in a perturbation expansion in ^{z) and a gradient expansion. 



Use is now made of the results from Section II B. Let A(<I>) be the energy of the lowest eigenstate of the transfer 
operator Htr + s^* + s*^ where Htr is given by (|19|) . The free energy per unit length of a single chain in the presence 
of a slowly varying field $(2) is 
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F,mz)] = yoA(O) - xi(0) (l * P -e? I ^ I') - X3 I $ I' (46) 

where Xi(0) is the static hnear susceptibihty of the system given by ( p^ ) and X3 is the third-order nonhnear suscep- 
tibihty given by (|3^). 

Near the three-dimensional transition temperature the continuum hmit perpendicular to the chains ($1(2) 
$(a:,y,z)) can be takenE2l. Assume that the interchain interaction Jy is non-zero only for nearest neighbour chains 
and that its value is Jx/4: and in the x and y directions respectively. Then the interchain interaction term in 
( p^ ) becomes 



^ J^^y dzRc[<i>,(z)*$^(z)] = -^-1- y dxdy M $ 



+ ^ 1^1 (47) 



where Ux and Cj, are the lattice constants perpendicular to the chains and J = ■^{Jx + Jy)- The final free energy 
functional ( p| ) is of the form (||) with coefficients 

^ = j-Xi(0) (48) 
B - -X3 (49) 
a = 1^ (50) 

= ^ (51) 

=xi(0)e?. (52) 
The three-dimensional mean-field temperature T^d is defined by the vanishing of the coefficient A, 

l = Jxi(0) (53) 
where the right hand side depends on temperature. This equation can be written in the dimensionless form 

J 1 



\a\ Af{K{T3D))' 



(54) 



The relationship between the interchain coupling J and the value of k at the three-dimensional transition defined by 
this equation is shown in Figure 0. 

IV. MEASURABLE QUANTITIES AT THE THREE-DIMENSIONAL TRANSITION 

It will now be assumed that the coefficients a, b, and c are temperature independent and T^d <^ Tq, so \a\ ~ a'. 
Consequently, the only temperature dependence in the fiuctuation parameter k is the factor of T- (compare equation 
(^). Scalapino, Imry and PincusB, MannevilleQ, and DieterichQ also made this assumption. It is an open question as 
to how realistic this assumption is for different microscopic models of the CDW transition. i_. 

First, this assumption allows us to rewrite n in terms of the width of the one-dimensional critical region AtiDt3 and 
evaluate the three-dimensional transition temperature as a function of the interchain coupling (Figure Roughly 

noV 1.5 J 



To J \a\Atljj- 
This is a useful relation because it gives a criterion 

J<^\a\Atljj (56) 
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for how weak the interchain couphng must be in order for the three-dimensional transition temperature to be sub- 
stantially less than the single-chain mean- field temperature Tq. 

The coefficient A' (which determines the specific heat jump and the coherence lengths) depends on the k dependence 
of /(k) and the T dependence of the fluctuation parameter k 



A' = 



1 dlnxi(O) 



J d\nT 
din \a\ 



c?ln/(K) dluK 



d\nT 



d In K d In T 



T=T3D 

Assuming that the coefficients a, b, and c are temperture independent ( [5^ ) reduces to 

^,^_2d\nf{K) 



J din K 



(57) 



(58) 



The specific heat jump, coherence lengths, and critical region width can now be calculated, ft turns out that they 
can each be written in a dimensionless form that only depends on the value of the fluctuation parameter k at the 
three-dimensional transition. The specific heat jump per chain given by (|^), ( ^ ) and (|5^) is 



AC: 



dn 



(59) 



where is the single chain coherence length given by (|f^). Hence the specific heat jump only depends on two 
quantities, smd k (or equivalently the interchain coupling) (see Figure |). Note that as the coherence length 
decreases the specific heat jump increases. For moderate interchain coupling (J > 0.05|a|) the specific heat jump 
becomes independent of the interchain coupling. The origin of this independence is not clear. 
The longitudinal coherence length is 



Co 



^0,1(1^) - Ao,o(k) 



dn 



-1/2 



(60) 



The dependence of the longitudinal coherence length on the interchain coupling is shown in Figure 

Quantitative comparison of these expressions for the specific heat jump and longitudinal coherence length with 
experimental data is not possible without accurate values of ^0 from microscopic theory. However, for all interchain 
couphng strengths the product of the specific heat jump and the longitudinal coherence length, AC^oz, depends only 
roughly on the interchain coupling (see the dashed curve in Figure 0) . Since this quantity is dimensionless it can be 
compared with experiment. If AC3D is the specific heat jump per unit volume then 



ACaDflajay^oz 



1.2-2.6. 



(61) 



For the blue bronze K0.3M0O3 experiments gives AC3d/A:b = (2.6±0.3) x lO^^l"^, $0z = 15±3 jM, and a^Uy = lel^ 



which gives AC3DCixa.yS,oz / ks = 0.6 ± 0.2. 



The transverse coherence length is 



£,0x 
CLx 



1/2 
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d\nf{K) 
din K 



-1/2 



(62) 



The dependence of the right hand side on k, is fairly weak. Consequently the transverse coherence length depends only 
weakly on the interchain coupling (Figure ||). This result is somewhat counterintuitive: it might be expected that 
as the interchain coupling increases the transverse coherence length increases. However, is not just determined 
by the coefficient Cx but also by the coefficient A' (^oa; = {Cx/A'Y/'^). It turns out that both these quantities have 
roughly the same dependence on J. PP 
A crossover temperature, Tx, can be defined at which correlations between chains become weak, by S,x{Tx) — axuxl. 
Hence, 



Tx 

T3D 



(63) 
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and it can be seen from Figure ^ that for a tetragonal crystal /Six ~ 0.1 virtually independent of any parameters. 
This value is comparable to the value obtained for by SchulzD for a slightly different model involving only fluctuations 
in the phase of the |Order parameter. This narrow crossover region is consistent with estimates for K0.3M0O3 based 
on X-ray scatteringlj. It may be possible to argue that this narrow crossover region, which is comparable to the 
width of the critical region estimated below, is consistent with the neglect of interchain fluctuations in the derivation 
in Section III. 

The width of the critical region given by the Ginzburg criterion (^) is 

. , _ ^ 9{i^?{x^A^) - XoA^)? .... 

Note that for a tetragonal crystal {Jx = Jy = J) this depends solely on the fluctuation parameter k. Figure [s] shows 
the resulting dependence of At^D on the interchain coupling. It is striking that the width of the critical region is 
independent of any parameters, Atsc ^ 0.05 — 0.08. As for the transverse coherence length this result is somewhat 
surprising since intuitively it might be expected that as the the interchain coupling increases the system becomes less 
one dimensional and fluctuations decrease. However, this independence is result of the transverse coherence lengths 
being only weakly depjendent on the interchain coupling. The values obtained for At^D are comparable to the value 
obtained for by Schulzl] for a slightly different model. It is interesting that even if the width of the one-dimensional 
critical region is large, say of order Tq (i.e., ATio ^ Tg), the actual observable critical region will be much smaller 
(AT3U <C T^r) <C To). This small width of the critical region is important because it shows that a mean-field treatment 
of the three-dimensional Ginzburg- Landau functional (||) is justified except in a narrow temperature region close to 
T^D- Note that this result is based on the assumption that the temperature dependence of the single chain parameters 
a, 6, and c can be neglected. i— . 

The value of At^D ~ 0.05 — 0.08 can be compared to experimental results. The analysis in ReferenceEJ and the 
data of Referencell3 suggests that At^^jj ^ 0.01 for the CDW transtion in K0.3M0O3. The data of ReferencalJ suggests 
that At^D ~ 0.01 for the SDW ttanstion in (TMTSF)2PF6. It should be stressed that these are very rough estimates 
based on the Ginzburg criterionc3. 

Rough estimate of n in several CDW materials. Since k is such an important parameter it is worthwhile making 
a very rough estimate of its value in various materials. It turns out that most materials lie outside the regime of 
validity of the lowest-level approximation, discussed in the Appendix. 



The longitudi nal coherence length given by (A15) can be used to estimate k. Equation (E3|) and the lowest-level 



approximation, (A7), can be used to provide a rough estimate of ^1 



1/2 

a =^ Co ( - ) (65) 



where ^0 is the single chain coherence length given by (|l5|). Hence, k is roughly given by 



(66) 

^20 can be estimated from X-ray and neutron scattering experiments (Table ||) . Deriving an accurate value of ^0 from 
microscopic theory is a subtle matter which will got be discussed here. For purposes of this paper it is sufficient to 
use the results of the simplest microscopic theoryiO. In that case 

^° " ^(oT ^''^ 

where Vp is the Fermi velocity and A(0) is the zero-temperature energy gap. Table || lists estimates of ^^o; 
and A(0) for several materials that undergo a three-dimensional CDW ordering transition. These values are used in 
equation (|6^) to give an order of magnitude estimate for k in each material. In all the materials k > 10^^ suggesting 
that the lowest-level approximation will give poor results. Given the experimental uncertainty in the parameters 
discussed above it is hoped that this discussion will generate more precise studies of this question. 

V. CONCLUSIONS 

The main points of this paper are the following. (1) To derive a reliable Ginzburg-Landau functional to describe 
three-dimensional ordering transitions in quasi-one-dimensional materials one must give a careful treatment of the 
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large intrachain order-parameter fluctuations. (2) Most physical properties are determined by a single dimensionless 
parameter k, defined in equation (^), which is a measure of the size of the intrachain fluctuations. (3) Commonly 
studied charge-density-wave materials lie outside the regime of the lowest-level approximation and so their accurate 
description requires an exact treatment of the intrachain fluctuations. (4) If the single chain Ginzburg-Landau 
coefficients are assumed to be temperature independent then the width of the three-dimensional critical region, 
estimated by the Ginzburg criterion, is virtually independent of any parameters, At^D ^ 0.05 — 0.08. 

This paper leaves a number of open questions and opportunities for future work, (a) It needs to be established 
whether the temperature dependence of the single chain Ginzburg-Landau coefficients a, b, and c is important in 
realistic microscopic models. If so how are the results of Figures |l|, ^ and H modified? (b) Quantitative comparison of 
the results here with measurable quantities such as the specific heat jump and longitudinal coherence length requires 
values of the single chain coherence length from microscopic theory, (c) For moderate interchain coupling the 
specific heat jump is given by AC^o — 0.7, independent of the strength of the interchain coupling. Is there a simple 
physical argument that can justify this simple result? (d) This paper treats the intrachain fluctuations exactly and 
the interchain interactions in the mean-field approximation. This neglect of the interchain fluctuations is presumably 
reasonable if the system is sufficiently anisotropic. It is desirable to have a quantitative criterion that defines the 
validity of this treatment. Furthermore, presumably when the interchain interactions become sufficiently strong the 
system is better described by an anisotropic XY model in which the fluctuations in all three directions are treated 
on an equal footing. It would be nice to have a description of the crossover from a transition at T^d ^ To (ensured 
by the criterion (pq)) to the case T^jj ^ Tq for strong interchain coupling. It is not clear how to attack this problem 
within any particular approximation scheme. It may be necessary to use a numerically intensive technique such as 
series expansions or a Monte Carlo simulation. Hopefully this paper will stimulate such a study. 



ACKNOWLEDGMENTS 



I have benefitted from numerous discussions with J. W. Wilkins. This work was stimulated by discussions with J. W. 
Brill. I am grateful to him for showing me his group's data prior to publication. I thank J. W. Brill, G. Mozurkewich 
and J. W. Wilkins for helpful comments on the manuscript. Some computational assistance was provided by M. C. B. 
Ashley. Some of this work was performed at The Ohio State University where it was supported by the U.S. Department 
of Energy, Basic Energy Sciences, Division of Materials Science and the OSU Center for Materials Research. Work 
at UNSW was supported by a Queen Elizabeth II Fellowship from the Australian Research Council. 



* electronic address: ross@newt.phys.unsw.edu.au 
^ G. A. Toombs, Phys. Rep. 40, 181 (1978). 

^ R A. Lee, T. M. Rice, and R W. Anderson, Phys. Rev. Lett. 31, 462 (1973). 

^ D. J. Scalapino, Y. Imry, and P. Pincus, Phys. Rev. B 11, 2042 (1975). 

*P. Manneville, J. Phys. (France) 36, 701 (1975). 

^ A. R. Bishop and J. A. KrumhansI, Phys. Rev. B 12, 2824 (1975). 

® W. Dieterich, Adv. Phys. 25, 615 (1976). 

H. J. Schulz, in Low-Dimensional Conductors and Superconductors, edited by D. Jerome and L.G. Caron (Plenum, New 

York, 1986), p. 95. 

* Proceedings of the International Conference on Science and Technology of Synthetic Metals, Goteberg, Sweden, August 

12-18, 1992, Synth. Met. 55-57 (1993). 
^ G. Griiner, Rev. Mod. Phys. 66, 1 (1994). 

^° F. Sharifi, A. V. Herzog, and R. C. Dynes, Phys. Rev. Lett. 71, 428 (1993). 

" M. Hase, L Terasaki, and K. Uchinokura, Phys. Rev. Lett. 70, 3651 (1993); J. P. Pouget, L. P. Regnault, M. Ain, B. Hennion, 
J. P. Renard, P. Veillet, G. Dhalenne, and A. Revcolevschi, Phys. Rev. Lett. 72, 4037 (1994); and references therein. 
See for example, S. Girault, A. H. Moudden, and J. P. Pouget, Phys. Rev. 39, 4430 (1989). [X-ray scattering]; J. P. Pouget, 
B. Hennion, C. Escribe-Filippini, and M. Sato, Phys. Rev. 43, 8421 (1991) [neutron scattering]. 
J. W. Brill, M. Chung, Y.-K. Kuo, X. Zhan, E. Figueroa and G. Mozurkewich, submitted to Phys. Rev. Lett. 
For accurate measurements of the specific heat near the SDW transition in (TMTSF)2PF6, see J. Coroneus, B. Alavi, and 
S. E. Brown, Phys. Rev. Lett. 70, 2332 (1993). 

L. D. Landau and E. M. Lifshitz, Statistical Physics, 2nd. ed., (Pergamon, Oxford, 1969), p. 478. 
^'^ D. J. Scalapino, M. Sears and R. A. Ferrell, Phys. Rev. B 6, 3409 (1972). 



11 



"R. H. McKenzie and J. W. Wilkins, Phys. Rev. Lett. 69, 1085 (1992). 

1* K. Kim, R. H. McKenzie, and J. W. Wilkins, Phys. Rev. Lett. 71, 4015 (1993). 

19 
20 
21 
22 



32 



R. H. McKenzie, Bull. Am. Phys. Soc. 39, 109 (1994). 

H. Castella, D. Baeriswyl, and K. Maki, J. Phys. (France) IV 3, 151 (1993). 
H. J. Schulz and C. Bourbonnais, Phys. Rev. B 27, 5856 (1983). 

The derivation of these coefhoioHtfi-js a subtle matter. For exaipele, for a CDW transition the presence of a pseudogap in 
the electronic density of stateilcSH and thermal lattice motioncJ significantly modify these coefficients. 
V. L. Ginzburg, Fiz. Tverd. Tela 2, 2031 (1960) [Sov. Phys. Solid State 2, 1824 (I960)]; S. K. Ma, Modern Theory of 
Critical Phenomena, (Benjamin/Cummings, Reading, 1976), p. 94. The width of the critical region, AT, is defined by the 
temperature at which the fluctuation contribution to the specific heat below the transition temperature, calculated in the 
Gaussian approximation, equals the mean-field specific heat jump AC. It should be stressed that this gives only a very rough 
estimate of the importanca-af fluctuations and that there are several alternative definitions of the width of the critical region 
(see for example ReferenceEll) . Consequently, care should be taken when comparing estimates from different references. This 
is particulary true since definitions in different references often differ by numerical factors such as 327r^! 
Z. Y. Chen, Phys. Rev. B 41, 9516 (1990). 

If a is given by ( |ll| ) then li-^ 2t^(Ati_D/(l — t))^ where Atio is given by the one-dimensional Ginzburg criterion (^ij). 
Scalapino, Sears and FerrelO assumed that T ~ Tq and in their notation At = 2Ati_D. Hence (1 — t)/At in their paper 
corresponds to (1/4k)^''^. 

The code used was an adapted version of Example 3 in S. E. Koonin and D. C. Meredith, Computational Physics, Fortran 
Version, (Addison Wesley, Redwood City, 1990), p. 67ff. 

See for example, E. Merzbacher, Quantum Mechanics, 2nd. Edition, (Wiley, New York, 1970), p. 457. 
J. M. Ziman, Elements of Advanced Quantum Theory, (Cambridge, Cambridge, 1969), p. 53ff. 
J. W. Negele and H. Orland, Quantum Many-Particle Systertts^ (Addison Wesley, Redwood City, 1988), p. 197ff. 
For more details on taking the continuum limit compare Ref.cJ, p. 208. 

M. R. Hauser, B. B. Plapp, G. Mozurkewich, Phys. Rev. B 43, 8105 (1991). In the notation of this reference (compare their 
equation (1)) the Ginzburg criterion gives Ar3D — 2{G/H)'^ ~ 2 K. It should be stressed that the width of the critical 
region, AT = 8 K, identified in this reference is based on an alternative criterion to the Ginzburg criterion^. 
M. J. Rice and__S. Strassler, Solid State Commun. 13, 1389 (1973). 
M. V. Sadovskii, Zh. Eksp. Teor. Fiz. 66, 1720 (1974) [Sov. Phys. JETP 39, 845 (1974)]; Fiz. Tverd. Tela 16, 2504 (1974) 
[Sov. Phys. Solid State 16, 1632 (1975)]. 

W. E. Masker, S. Marcelja, and R. D. Parks, Phys. Rev. 188, 745 (1969). 
J. R. Tucker and B. I. Halperin, Phys. Rev. B 3, 3768 (1970). 



FIG. 1. Dependence of the three-dimensional transition temperature Tso and the fluctuation parameter k on the interchain 
coupling J. It has been assumed that the single chain Ginzburg-Landau parameters a, b, and c are independent of temperature. 
The reduction of Tso below the mean-field transition temperature Tb for a single chain also depends on the width of the 
one- dimensional fiuctuation region Atio- 



FIG. 2. Dependence of the specific jump per chain AC and the longitudinal coherence length ^oz on the interchain coupling 
J. It has been assumed that the single chain Ginzburg-Landau parameters a, b, and c are independent of temperature. The 
longitudinal coherence length is normalized to the single chain coherence length ^o defined by (^5|). Note that for moderate 
interchain coupling (J > 0.05|a|) the specific heat jump is determined solely by the the single chain coherence length ^q. 



FIG. 3. Very weak dependence of the width of the critical region for the three dimensional transition and the transvepe 
coherence length on the interchain coupling. The width of the critical region At^n is calculated from the Ginzburg criterion^. 
The transverse coherence length ^ox is normalized to the lattice constant ax- The results shown are for a tetragonal crystal. 
For an orthorhombic crystal the vertical scale is changed by a factor depending on the anisotropy (see equations (^2|) and (^)). 
This figure establishes the striking result that the width of the critical region is virtually parameter independent. 



FIG. 4. Failure of different approximation schemes for treating the effects of fiuctutations in a single chain. The dependence 
of the order parameter susceptibility function /(ac) (defined in equation (^^) on the fiuctuation parameter k (defined in equation 
(^)) is shown. The solid curve is the exact result obtained by numerioally solving the eigenvalue equation (^8|). The datted 
curve is the result, /(k) = 1/k, of the lowest-level approximation (LLA)u which has been-nsed in previous tpatments ofQ'u'Q of 
quasi-one-dimensional phase transitions. The dashed curves are the results of the HartreetJ and Hartree-FoclO approximations. 
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FIG. 5. Deviation of the lowest-level approximation from the exact results. The linear susceptibility /(k), the third-order 
susceptibility <;(«;) (defined in equation (^H])), and the derivative ^}"^ (which determines the coefficient A' (see equation (^)), 
all normalized to their values in the lowest-level approximation, are shown as a function of the fluctuation parameter k. Note 
that the horizontal scale is logarithmnic. Clearly the lowest-level approximation is only quantitatively reliable for k < 10"''. 
In Table ^ it is estimated that typical materials have 10~^ < k < 1. In this range the lowest-level approximation will predict 
a specific heat jump up to an order of magnitude smaller than the exact results. 



TABLE I. Summary of symbols for the important quantities for a single chain and a three-dimensional (3D) system of 



weakly coupled chains. The items 


in parentheses refer to relevant equations and figures. 




Quantity 


Single chain 


3D 


Order parameter 

Ginzburg-Landau 
free energy 
Meau-field 

transition temperature 

Ginzburg-Landau 

coefficients 

Specific 

heat jump 

Intrachain 
coherence length 

Fluctuation 
region width 


t 

To 

a, b, c 
ACiD ^ 

^0 (0) 
AtiD (0) 


$(a;,y (g) 
# 

nn < To [Fig. |l[ 
A, B, C.,Cy, a 
AC [(H), Fig. 1 
Co. [(0), Fig. 1 
At3D [(g), Fig. 1 



TABLE II. Order of magnitude estimate of the fluctuation parameter k in several materials that undergo a three-dimensional 
charge-density-wave transition. The intrachain coherence length ^oz, measured by X-ray or neutron scattering, and the sin- 
gle-chain coherence length are used to estimate the value of the fluctuation parameter k at the three-dimensional ordering 
transition (see equation (|6^)). is estimated from the Fermi velocity vp and the zero-temperature energy gap A(0), using 
equation (p7t), based on the simplest microscopic theory. The estimates of foz, vf, and are probably only accurate to within 
a factor of two. Consequently, the estimate of k is only accurate within an order of magnitude. Nevertheless, clea rly a ll these 



A 2 ) and so 



materials lie outside the regime of validity (k < 10~ , see Figure q) of the lowest-level approximation (see Section 
an exact treatment of the intrachain fluctuations of the order parameter is necessary. Also note that these values of k can be 
used to estimate the values of J/|a| and they are in roughly the same regime as considered in Figures 3-5. 



Material 




10 m/sec 


A(0) 
meV 




K 


K0.3M0O3 


IS'' 


2.0*^ 


100^= 


3 


0.2 


(TaSe4)2l 


60'' 




200' 


10 


0.1 


K2Pt(CN)4Bro.3 


100^ 


10'' 


100' 


15 


0.1 


TSeF-TCNQ 


30j 


0.15^ 




2 


0.02 
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APPENDIX: PREVIOUS APPROXIMATE TREATMENTS OF FLUCTUATIONS 



1. Hartree and Hartree-Fock treatment 

For completeness these approximations which have been considered in the Uterature0@'0 are discussed here. The 
replacement | (f>{z) 2q <\ (j) p>| (/'(z) P is made in_tlie free energy functional (|l|) in the partition function 
(IitI). The cases q = 1 and q = 2 correspond to the Hartrea^ and Hartree- FockEj approximations, respectively. The 
functional integral is then over a Gaussian field and can be performed analytically. In effect one makes the Gaussian 
approximation with 

a ^ a + T, = a + 2qb <\ (f>\^> . (Al) 
The expectation value <| 4> P> is calculated self consistently 

=<| ^ |2>^ _ / dk = , (A2) 

2qh 2^ J a + E + cfc2 2Vc(a + S) 

Consequently, the dimcnsionless self energy E = S/|a| satisfies the self-consistent equation 

^-i^S'^ (A3) 



Since the linear susceptibility is 



2(1] -1) 



xm-^ (A4) 



the function /(k) in the susceptibility is given by 

/(k) = -J^ . (A5) 

The results of the Hartree and Hartree-Fock approximations are shown as dashed curves in Figure H; they are 
in qualitative but not quantitative agreement with the exact results. As k — * 0, S ^ 1, and so ( |A3[ ) implies 
S ^ 1 + q^Kl2. Hence, this approximation gives 

/W-^TT^ as K^Q. (A6) 



1. The lowest- level approximation 



Close to the three-dimensional ordering transition k ^ 1. This limit corresponds to the semi-classical limit of the 
Schrodinger-type equation (^8|). Previously a Gaussian wavefunction, sharply peaked at p = 1, has been used as a 
variational ground state wave function for (|2^). The results areo 

\yi - Ao,o = (A7) 
In this approximation the correlation length (given by (||), (|2^) and (|A7|)) is 

a(T) = ^. (A8) 

This result was given previously in References!, 0, and^. In this limit it is to be expected that has an expectation 
value of 1 in the ground state. If it is further assumed that the n = states form a complete set (the lowest-level 
approximation) then 



1 =< uo^oI/O^Imo.o >=< uo,a\pe } |uo/ >< uqj\ pe*'^|wo,o >=|< woa|p|uo,o >|^ (A9) 
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The expression ( |3l| ) then reduces to 



/(«) = z- 



(AlO) 



Note that in the lowest-level approximation the sum on the left hand side of (|3^) is only 1/2. 

In the lowest-level approximation the dimensionless function g{n) in the third-order susceptibility is given by 



(Ao,i - Ao,o)^ 



(Ao,i — Ao,o)^(Ao,2 — Ao,i) 



7 

2^ 



(All) 



Figure H shows how thegesults of ( |A10|) and ( All ) deviate significantly from the exact results for k > 10^'^. Previously 
Bishop and KrumhanslQ pointed out the shortcomings of the lowest-level approximation for the case of a real order 
parameter. _ _ _ 

The three-dimensional transition. Scalapino, Imry and Pincug3, Mannevilleu, an d DieterichO studied the three 
dimensional transition using the lowest-level approximation (defined in Section A 2), i.e., they assumed that only 
the lowest eigenstate (for n = 0) for each of the angular momentum values £ = 0, ±1, and ± 2, is important. As 
mentioned earlier they also assumed that the single chain coefficients a, 6, and c are independent of temperature. The 
zero-wavevector linear susceptibility is given by 



XI (0) 



Solving (p3|) for the transition temperature gives 



2a^c 



At I a I 52T2 ■ 



(A12) 



T^D — ■v 

Equation ( ^9| ) for the specific heat jump reduces to 

AC 

The coherence length parallel to the chains is 

The transverse coherence lengths (for = Jy) are 

^ _ ax 



'2Jc 
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(A13) 



(A14) 



(A15) 



(A16) 



The above expressions ( |A13 - A16) were previously obtained by Mannevillea and DieterichO. 

In the lowest-level approximation the expression (^4|) for the width of the critical region (for = Jy) reduces to 



0.08. 



(A17) 
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